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Abstract. 

CZJ , Recently, we predicted based on symmetry considerations that KMnFeFg and 

Ba2Ni7Fi8 are likely magnetoelectric multiferroic materials. In this contribution, we 
investigate with Landau theory and crystal structure considerations the polarization 
and the linear magnetoelectric effect in these materials. Based on these two examples, 
we show that any magnetoferroelectric will display additional electrical polarization 
below its magnetic ordering temperature. This additional electrical polarization is not 
related to the linear magnetoelectric effect. Its magnitude depends on the dielectric 
'^ ' susceptibility. 
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1. Introduction 

In recent years, the coupling between magnetic and dielectric properties in transition 
metal oxides gave rise to a significant research effort [U l2l [3]. This effort is governed by 
the emergence of new fundamental physics and potential technological applications [21 13]. 
This recent research is mostly concentrated on oxides |lj. We think it is very important 
to search for new materials, especially for materials different than oxides [5l [6] . Since 
multiferroic and magnetoelectric properties rely on symmetry considerations, the use of 
group theory is of prime importance to look for new materials. The use of symmetry 
arguments in order to predict new multiferroic or magnetoelectric materials is not new 

In this contribution, we use the phenomenological Landau theory to describe 
the previously predicted magnetoelectric effect in multiferroic fluorides KMnFeFg and 
Ba2Ni7Fi8 [5]. These two materials exhibit a pyroelectric symmetry and thus exhibit 
a spontaneous electrical polarization at room temperature. Here we are interested in 
the behavior of this electrical polarization at the magnetic ordering temperature. For 
this reason, we will consider in the free energy expression only the lowest degree term 
of the electrical polarization which is P^. The other coupling terms between the various 
order parameters are derived from symmetry considerations. Every product of the order 
parameters components belonging to the same irreducible representation is invariant by 
time reversal and the crystal symmetry |10j . 

We use the general case of Ba2Ni7Fi8 (symmetry PI) in order to discuss various 
contributions to the induced polarization. In particular we show that any magnetically 
ordered pyroelectric materials will display additional electrical polarization below its 
magnetic ordering temperature. This electric polarization is not related to the linear 
magnetoelectric effect. We discuss these results in the light of recent reports on 
magnetoelectric and multiferroic materials [HI [121 (El dH fT5] . 

2. Study of KMnFeFe 

The crystal structure of the fluoride KMnFeFg has been discussed in detail in literature 
[T6] . This compound crystallizes in the space group Pba2 (n°32), where the Mn and Fe 
ions order on the 4c Wyckoff position of the structure and occupy statistically the 4b 
Wyckoff position. KMnFeFg orders ferrimagnetically below T^^ = 148 K with a ratio 
^=3. The magnetic unit-cell is identical to the chemical unit cell and thus A; = . 
A symmetry analysis by Bertaut's method gives rise to the results presented in Tables 
[I]and[2][l6l[n!. 

For the Wyckoff position 4c, the labelling of the various spins is Si in (x,y,z); 5*2 
in (x,y,z); S3 in (l/2-x,l/2+y,z) and ^4 in (l/2+x,l/2-y,z). The authors defined the 
following magnetic vectors: 
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Magnetic space groups 


Ti 


G. 


Ay 


c. 


Pba2 


Ts 


c. 


Fy 


G, 


Pha'2' 


Ts 


A. 


G, 


F, 


Pb'a'2 


r4 


F. 


Cy 


A, 


Pb'a2' 



Table 1. Irreducible representations for the space group Pba21' associated to k=0 for 
the Wyckoff position 4c [T6j . 
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Table 2. Irreducible representations for the space group Pba21' associated to k=0 for 
the Wyckoff position 2b [TB] . 



F = Si -\- S2 -\- S^ -\- S4 
G ^ oi — 'J2 "I" '-'3 — '-'4 

G = 61 + 62 — 03 — 64 
A =61 — 02 — 03 + 64 



(1) 



For the Wyckoff position 2b, the authors defined two magnetic vectors F' = 5'i + 5'3 
and C = Si — Ss where this time Si is in (0,l/2,z) while 5*3 is in (1/2, 0,z). The 
decomposition is given in Tables [1] and [2l The neutron data show that the best model for 
the magnetic structure is given by the r4 mode. The various components of the magnetic 
vectors defined in ([1]) transform as the reducible representation. The appropriate way 
is to consider vectors which transform as irreducible representation (IR) [T^. In order 
to simplify the expression of the free energy and fulfill the symmetry requirements, 
we defined the magnetic order parameter rj which transform as the r4 IR as observed 
experimentally. The total magnetization of the system is defined by the order parameter 
M. Its components transform like the components of F (tables [1] and [2]). Thus M^^, M^ 
and M^ transform like the IRs r4, r2 and T^ respectively. 

7] and M^ transform both as the r4 IR, thus the lowest coupling term between r) 
and M IS rj M.^- The crystal structure of KMnFeFg has an orthorhombic pyroelectric 
symmetry. Thus this material exhibits a spontaneous electrical polarization along 
the z axis. All the symmetry elements of the space group Pba2 let invariant only 
the P^ component of the electrical polarization. Thus the lowest coupling term 
between the magnetic order parameter r] and the polarization is rj^Vz- For the possible 
magnetoelectric coupling terms, one has to consider the possible terms of the type rjMP 
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j. It can be shown that the electrical polarization components P^;, P^ and P^ belong 
respectively to the IRs Ps, P4 and Pi respectively. Looking at the possible rjMP and 
taking into account the above considerations, we can write the free energy for KMnFeFe 
in the presence of a magnetic field as: 



+ Y, ^Mf + XivPyM, + X2vPzM^ - M.H 
We find the partial derivatives of $ at the equilibrium conditions. 
(9$ 



(2) 



(9$ , 

— = ar] + br]' + 7M, + 2ar]P, + XiM^Py + AsM.P, = 

= KnP. = 

dPx 

-^ = f^22Py + Air^M, = 

— - = KgaP. + \2r]M. + ar/' = (3) 

oPz 

<9$ 






= CuM^ + X2r]Pz + 'j'n- H^ = 







— = c,.,M, -H, = 



9$ 

Ml 



C33M, + XiriPy -H, = 



In order to investigate the linear magnetoelectric coupling, we express the electrical 
polarization as function of the magnetic field. For this purpose, we use the expressions 
of ^ and of ^. We find: 

Py = ^vMz 

1^22 (^^ 

_ -{ay^ + X2VM,) ^ ^ 

-'2 

Using the results of equations ([3]) and (jlj), we can derive the expression for the 
various components of the magnetization. 

,. ^^33 rr , v{^2(rrf - l^-,3l) 

%3Cll - X^^T %3Cll - X^2T 

My = ^ (5) 

C22 

M. = ^^^V^^. 

K22C22 - Xfi]^ 
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From equations (jl]) and i^, we can find the expression for tlie linear relationsliip 
between tlie induced polarization and tfie application of a magnetic field (linear 
magnetoelectric effect) for the fluoride KMnFeFg. 



' K22Css-Xiv' (6) 

_ -A2?? (A27 - (TCuW 

-'2 ^ o rt-tlx ~r 



fi^sscii - A^?]^ ^ K33C11 - Xli]'^ 

From equation (Q, we find that there are two non-zero components for the linear 
magnetoelectric tensor of KMnFeFg: 023 and a^i. 023 takes the value ~ ^'^a 2 and 
ctsi is equal to ~ ^^^ 2 • We notice that the polarization along x is zero at any 
temperature even under the application of a magnetic field. Our results suggest that 
the polarization along y is purely induced by the magnetic field (linear magnetoelectric 
effect) while the polarization along z would have an additional component arising at 
the magnetic transition temperature (term in rj'^ in equation ([6])). We notice that this 
additional electrical polarization along z arising below Tjv is present even in the absence 
of a linear magnetoelectric effect (Aj = 0). We stress that our model describes only 
what is happening below the magnetic ordering temperature. One should keep in mind 
that there is already a spontaneous polarization along z at room temperature resulting 
from the pyroelectric symmetry exhibited by KMnFeFg. We notice that the magnetic 
point group m'm2' described by the irreducible representation r4 should exhibit two 
non-zero components for the magnetoelectric tensor according to Ref. [T9j: a23 and 
^32 instead of 031. This difference can be explained by the fact that we have used the 
labelling of the authors of Ref. [16]. If we inverse S3 and 5*4, we find that the total 
magnetization component M^^ transforms like r4 instead of r2 (see tables [1] and [2]). 
Consequently the invariant term responsible for the linear magnetoelectric effect is not 
\2rjPzMx but \2r]PzMy which will give rise to 032 instead of 031. In addition to this 
linear magnetoelectric effect, we notice that a spontaneous magnetization component 
along X can be displayed below the Neel temperature in the absence of magnetic field 
(see Eq. (ED). 

3. Study of Ba2Ni7Fi8 

3.1. Estimation of the spontaneous polarization^ 

To estimate the spontaneous polarization of a pyroelectric system, one may deduce 
the possible high-symmetry (or high-temperature) structure from the low-symmetry 

§ We notice that this estimation is not possible for KMnFeFg since in the reported crystal structure 
[16j . the authors fixed the cations on the atomic positions determined in the PA2bc symmetry and not 
in Pba2. Consequently a redetermination of the crystal structure of KMnFeFg is required to investigate 
properly the coupling between dielectric and magnetic properties. 
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structure based on a pseudosymmetry analysis. This concept gave rise to the prediction 
of a large number of displacive ferroelectrics [20l [21] . It has been implemented in the 
program PSEUDO (Crystallographic Bilbao Server) [21]. When the atomic coordinates 
of a given structure display an approximate symmetry in addition to the actual space- 
group symmetry, the structure can be considered as pseudosymmetric with respect to 
a supergroup containing this additional symmetry. The existence of pseudosymmetry 
in a crystal structure indicates a slightly distorted structure of higher symmetry. If the 
distortion is small enough, one can expect the crystal to acquire this higher symmetry at 
a higher temperature. Using the crystal structure of Ba2Ni7Fi8 reported by Lacorre et al. 
[22] , we performed a pseudosymmetry search among all the minimal supergroups of PI 
(symmetry exhibited by Ba2Ni7Fi8). Of the minimal supergroups of PI, we determined 
PI was the only pseudosymmetric minimal supergroup. We therefore concluded the 
high-temperature paraelectric phase of Ba2Ni7Fi8 is in this space group. However, this 
pseudosymmetric minimal supergroup would produce a maximum atomic displacement 
of 1.3 A. The high value of the atomic displacement towards PI suggests that the 
paraelectric phase will not be reached before melting or decomposition. However, the 
possibility of twinning and domain switching may exist. The pseudosymmetric minimal 
supergroup allows us to estimate the spontaneous polarization in Ba2Ni7Fi8. This is 
can be done using an ionic model by simplifying the electric charge of each ion with a 
point charge [23] . 

Ps = [^)Y.rniQAzi (7) 

i 

e is the elementary charge, V the unit-cell volume, nij the multiplicity of the 
ion, Qj the ionic charge and Azj the displacement ion along the polar axis. Having 
discussed already the hypothetical high temperature paraelectric phase, we can estimate 
the displacement Azj for each ion. In that case, since we have the symmetry PI, we 
have to deal with a spontaneous polarization vector having three non-zero components: 
Pj,., ^y and P^. Using the point charge model, we estimate that P^; = 0.246 ^C/cm^; 
Vy = 10.87 tiC/crn^ and P^ = -26.87 ^C/cm^. The total polarization is thus ||^|| = 29 
/iC/cni^. This is a large spontaneous polarization but this value is an upper bound. The 
ionic model used is a crude picture which does not take into account the covalency of 
the bonds. Additionally, the calculation is strongly dependent of the accuracy of the z 
coordinate. For BaAl204, the ionic model gives an overestimation by a factor of six of P^ 
compared to the experimental value [23] . Nevertheless, the spontaneous polarization in 
Ba2Ni7Fi8 remains high and comparable to known multferroic systems such as YMnOa 
(P, = 5.5 fiC/cm'^) ^. 



3.2. Landau phenomenological description of the linear magnetoelectric effect 

In the previous fluorides, the magnetic frustration appeared in corner- sharing octahedra 
through the existence of interconnected triangles of magnetic ions. This leads to a 
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single type of magnetic exchange interaction. P. Lacorre and coworkers have been also 
investigating compounds like Ba2Ni3Fio (n = 9) and Ba2Ni7Fi8 (n = 21) which are 
members of the Ba6M„Fi2+2n family [221 [26] . In this family where M = Ni, there 
are not only corner-sharing octahedra but also edge-sharing octahedra. Both types of 
interaction exist in the Ba2Ni3Fio and Ba2Ni7Fi8 compounds. These compounds have 
been investigated by means of powder neutron diffraction at room and low temperatures. 





hi 


K 


Ti 


1 


-1 



Table 3. Irreducible representation for the space group Pll' associated to k=(0, 0, 
0). 

In order to keep consistency in the notations between the two cases, we will call 
also 77 the magnetic order parameter describing the magnetic ordering of Ba2Ni7Fi8. 
Since there is only one IR in that case, following similar arguments that in the case of 
KMnFeFg; all the coupling terms r]Mj are invariant. Since the crystal structure is PI, 
spontaneous electrical polarization components are non-zero for all the three directions 
X, y and z. Consequently, the invariant terms involving the electrical polarization 
are ^. • ^P^ and 'Yli^iV^Pi- ^^ addition to these terms, one needs to consider the 
magnetoelectric coupling terms. All the electrical polarization components transform 
as the unique Fi IR. Thus all rjPiMj are allowed. Consequently we can write the free 
energy for this fluoride as: 



^ = *0 + 1??^ + ^V"^ + {(TxPx + CTyPy + O.P,) 7]"^ 






+ (7,M, + ^yMy + 7,M,) r/ + (AiP,.M, + X^P^My + A3P.M,) t] (8 

+ {XiPyM^ + X^PyMy + AePyM, + A7P.M,) r/ 
+ (AgP.M, + AgP.M,) 7] - M.H 

M and P designate respectively the total magnetization and the electrical 
polarization of the system. Calculating the derivatives of the different variables, we 
can write the following set of equations: 
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_ =ar] + br]^ + 2r] {a^P^ + OyPy + a,P,) + 7^M^ + -fyMy + 7,M, 

+ AiP,.M, + A2P..M^ + X^P^M, + A4Pj;M, + X^PyMy + AeP^^M, 
+ A7P.A4 + AgP.Mj, + AgP.M, = 

a,?]' + n^P^ + (AiM, + AaMy + A3M,) r; = 



__ = ayri^ + /tj^Pj, + (A4M, + AsMj, + AgM,) r/ = 
__ = a,r,^ + /€,F, + (A7M, + X^My + AgM,) 77 = 
^^ = c^M^ + 7^77 + (AiP, + AsPj, + A3P.) r]-H^ = 



(9) 



dMy ""'' 



My + 7j,r^ + (A4PX + AsFj^ + AgP,) v - Hy = 



CzM, + 7,77 + (A7P. + AgPy + AgP,) i]-H, = 



dM, 



Here the dielectric tensor Kj ,,• has 9 non zero terms. To simphfy the notation, 
we write i^x= 1^11+ 1^12+ i^i'ii i^y= 1^21+ 1^22+ 1^22, and ^^.z=l^?,l+ 1^32+ 1^2,2,- AVe use a similar 
notation for Cx-, Cy and c^.In a similar manner as in the case of KMnFeFg, we can 
extract from equation ([9]), the expression of Mj as function of the magnetic field, of 
the magnetic order parameter and of the polarization. We find the expressions for the 
various magnetic components: 



My = "" '"'" ' '"7 ' '"^" ' "°^^^ (10) 

Cz 

If we replace the various Mj components in the expressions of ^, we will find, for 
instance, P^; as function of Py and P^. To determine the expression of P^; as function 
only of the various constants and of the magnetic order parameter rj, we need to express 
Vy and P2 as function of P^^. After we can use these results to solve ^p- = as the only 
function of P^;. We have in total six equations with six variables. This system is solvable 
exactly. However, the full expression of the various electrical polarization components 
is very lengthy and complicated. Consequently, in order to simplify the discussion we 
made a series expansion of P^., Vy and P^ in terms of r] at the second order: 



Hx- 


-V{lx + AlP. + A4Py + A7P.) 


Hy - 


Cx 
-V{ly + A2PX + X5Py + AgP.) 


Hz- 


Cy 

-l]{lz + A3P. + XePy + XgPz) 
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Px — ( H;j. H Hy H H^ 

\ ^x '^x ^y '^x ^z '^x 

\ '^x C-x rix Cy Hx (^z '^x / 

Py = ( Hx H Hy H H^ 



y riy Cxi^y Cy'^y C^z^y J 

and P.= (^Hx + ^Hy + ^H.^ 

\C-xl^z CyKz CzKz 

V ^z CxKz CyKz CzHz / 

From the results of equation (TTTj) . we are able to determine the magnetoelectric 
tensor [oij]. All the terms of the tensor are non-zero. This is in agreement with the 
results of Ref. [19]. We find a similar result than in the case of KMnFeFg where the 
various component of the tensor are equal to —^^rj. This result is equivalent to the one 
reported in section [2] in the hypothesis of small value of r]. In a similar way, there are 
three non zero spontaneous magnetization components along the three directions x, y 
and z (Mj 7^0 for Hj 7^0 in Eq. (ITOi) ) To this induced electrical polarization induced by 
the magnetic field (linear magnetoelectric effect), there is another contribution of the 
induced electrical polarization which results from the coupling to the magnetic order 
parameter (terms in ^77^). We shall discuss these two induced electrical polarization 
contributions in section HI 

4. Discussion 

Using phenomenological Landau analysis, we have shown that for any magnetically 
ordered pyroelectric materials (= multiferroics) an additional electrical polarization 
should arise upon cooling trough the magnetic ordering temperature. In previous 
Landau theoretical treatment of magnetically ordered ferroelectrics [llj, terms linear 
in the electrical polarization and quadratic in the magnetic order parameter were not 
taken into account. However this kind of terms is allowed for any directions along 
which a spontaneous polarization exists. In various publications terms quadratic in 
the polarization and in the magnetic order parameter are considered [TTl IT2l [T3l fT¥[ flS] . 
These kind of terms are expected to have a lower contribution (due to their lower degree) 
compared to terms like ^^^Piff where rj is the magnetic order parameter, k the inverse 
dielectric susceptibility and Pj the direction of the spontaneous electrical polarization. 
We believe that this kind of term is responsible for the reported additional electrical 
polarization for YMnOs and LuMnOs below their respective Neel temperatures [271128] . 
Lee et al. explained their results by claiming that they observed the simultaneous 
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condensation of three order parameters, spin (S), lattice (L) displacement, and electric 
dipole (D) moment [27]. They claimed that -aS'^L'^ and -L.D terms condense at the 
same temperature, with T^ = T^ for a sufficient strong spin-lattice interaction strength 
a. If one considers the literature [HI [131 [IS] ^^^ our results, we see that this statement 
is incorrect. 

The results given by equations (jl]) and (ITTl) are similar. If we make a series 
expansion of P^^ and P^ in terms of r/ to second order, we find analogous expressions 
for both expressions. While the magnetic order parameter i] varies like a/T — Tjv (mean 
field approximation), the additional polarization resulting from the coupling term aiPiff 
varies linearly with the temperature (it is function of rf"^ see Eqs. (|6])and ( ITTl) ). More 
importantly, this additional polarization arising in the magnetically ordered phase is a 
function of the dielectric susceptibility. Consequently, the higher will be the dielectric 
susceptibility of the material under consideration, the higher will be the additional 
electrical polarization in the magnetically ordered phase (see equation (fTTj) ). This 
behavior differs from the linear magnetoelectric effect since it does not depend on the 
magnetic susceptibility [29]. In the light of these considerations, we claim that the 
observations made for YMnOs and LuMnOs are not isolated cases but should be true 
for any magnetically ordered pyroelectric materials. 

5. Conclusion 

In conclusion, we use a phenomenological Landau analysis to describe the linear 
magnetoelectric effect in two multiferroic fluorides, namely KMnFeFg and Ba2Ni7Fi8. 
We use the general expression of the free energy of Ba2Ni7Fi8 to discuss in detail the 
various contribution to the induced polarization. We show that contrary to previous 
reports the biggest contribution of the additional polarization below the magnetic 
ordering temperature arises from a term a^rfVi where r] is the magnetic order parameter 
and Pj the direction of spontaneous electrical polarization. This additional polarization 
arising below the magnetic ordering temperature is characteristic of any magnetically 
ordered pyroelectric material. We show that this induced electrical polarization is 
proportional to the dielectric susceptibility of the material. Thus one should see a 
change in the slope of the temperature behavior of the electrical polarization at the 
magnetic transition with an increase of its absolute value. We expect that this work 
will stimulate experimental investigations of the above reported fluorides but also of the 
crystal structure of the magnetoferroelectric materials below their magnetic ordering 
temperature. 
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